Collective Excitations of a Bose-Einstein Condensate in an Anharmonic Trap 
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We investigate the collective excitations of an one-dimensional Bose-Einstein condensate with 
repulsive interaction between atoms in a quadratic plus quartic trap. With using variational ap- 
proaches, the coupled equations of motions for the center-of-mass coordinate of the condensate and 
its width are derived. Then, two low-energy excitation modes are obtained analytically. The fre- 
quency shift induced by the anharmonic distortion, and the collapse and revival of the collective 
excitations originated from the nonlinear coupling between the two modes, are discussed. 
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I. INTRODUCTION 

One of the most important characters of an interact- 
ing quantum many-body system is its response to ex- 
ternal perturbations, where collective excitation modes 
represent a very effective tool for exploring the role of in- 
teractions and testing theoretical schemes. For this rea- 
son, measurements of the collective modes in the trapped 
gases of alkali-metal atoms 0, [E were carried out 
soon after the discovery of Bose-Einstein condensates 
(BECs). For the dilute degenerate gases, the essential 
physics of the BECs ground state is included in the 
Gross-Pitaevskii equation (CPE). The nonlincarity, orig- 
inated from the interatomic interaction, is included in 
the equation through a mean field term proportional to 
condensate density. The study about the collective exci- 
tations of the condensation has been investigated exten- 
sively by using various theoretical methods [J, IE IE 0> HI • 
The remarkable agreement between measured frequen- 
cies and theoretical predictions is one of most impor- 
tant achievements in the investigation of these new sys- 
tems. Meanwhile, because of the emergence of nonlin- 
earity, a lot of interesting phenomena in the collective 
excitations of BECs, such as frequency shift yj, mode 
coupling [E ^3, [nl llpl. damping 0, f3 > collapse and 
revival of oscillations |l5l fl6| and the onset of stochastic 
motions for strong driving amplitude [TtI fl8l. IT^ | , have 
been paid much attention. 

Recently, the collective dynamics of an one- 
dimensional trapped ultra-cold Bose gases has attracted 
considerable attention, since experiments on trapped 
Bose gases at low temperature have pointed out the oc- 
currence of characteristics ID feature [^E l2ll |2^ . In 
previous works the studies about the collective excita- 
tions of BECs in the magnetic traps are mainly limited 
to the harmonic case. However, in practical situation 
of experiments, the trap usually is not purely harmonic. 
With this concern, we study the collective excitations of 
the BECs in an one-dimensional in a harmonic trap with 
a quartic distortion. Our aim is to understand how the 



distortion affects the collective excitations of BECs. 

Our study is facilitated by variational approaches. Us- 
ing a Gaussian trial function, the GPE is transformed 
into a set of second order ordinary differential equations 
about some parameters that characterize the condensate 
wave function. Then we derived the expressions of the 
two low-energy oscillation modes analytically and the 
nonlinear coupling between the two modes is revealed. 
In particular, we find that, a very small anharmonic dis- 
tortion may cause a significant frequency shift of the 
excitation modes when the atomic interaction is strong. 
Finally, we demonstrate that the anharmonic distortion 
may give rise to the collapse and revival of the collective 
excitations. 

The paper is organized as follows. In Sec. II we derive 
the governing equations for the center-of-mass coordinate 
of the condensate and its width. In Sec. Ill we discuss 
the collective modes and the frequency shift caused by 
the anharmonic distortion. In Sec. IV, we demonstrate 
the collapse and revival of the collective excitations in 
anharmonic potential. Final section is our conclusion. 



II. VARIATIONAL APPROACH AND 
GOVERNING EQUATIONS 

We consider dilute degenerate bosons confined in a 
cigar shaped trap and assume that the system is far from 
the Tonks-Girardeau regime 23]. Then, the BECs can 
be well described by the dimensionless ID GPE, 



1 d 2 
2dx^ 



V(x)+g\ij(x,t)\ 2 }iP(x,t), (1) 



where the coordinate x is measured in unit of ^J%/mui x 
and time is in unit of 1/lu x . lu x is the x component 
frequency of the harmonic trap. ifj(x,t) is the macro- 
scopic wave function of the condensate normalized so 
that J \tp(x)\ 2 dx = 1; g = 47V Trails/ \/h/mu> x char- 
acterizes the interatomic interaction and is defined in 
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terms of the s-wave scattering length a s (below we shall 
be concerned with repulsive BECs for which a s > 0); 

aid = J \(p(y,z)\ 4 dydz/ (J \ip(y, z)\ 2 dydz) 5/2 is a coeffi- 
cient which compensates for the loss of two dimensions 
pij . In the above expressions, N is the total number 
of atoms and <p(y, z) is the ground wave function of the 
lateral dimensions. 

The trapping potential we consider takes form 



V(x) = \tf + \x*). 



(2) 



The quartic term in the potential denotes the anhar- 
monicity of the trap. In [2q. the authors created such a 
quartic confinement with a blue-detuned Gaussian laser 
directed along the axial direction. In their case, the non- 
rotating condensate was cigar shaped and the strength 
of the quartic admixture was A ~ 1CP 3 . Here, we re- 
gard A as a controllable parameter and assume that the 
anharmonicity is weak, i.e., |A| <C 1. 

The problem of solving Eq.© can be restated as a 
variational problem corresponding to the minimization 
of the action related to the Lagrangian density |f| 
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where the asterisk denotes a complex conjugate. In order 
to obtain the dynamics of the condensate in the trapping 
potential we will find the extremum of Eq. © with a set 
of trial functions. In our case, a natural choice of the 
trial function is a Gaussian, i.e., we take 



ip(x,t) = r}(t)e 
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The first one corresponds to the normalized condition of 
the wave function, J \4>(x, t)\ 2 dx = 1. Therefore, once 
we know the behavior of the center and width of the 
condensate, we can calculate the evolution of the rest 
of the parameters, and then completely characterize the 
dynamics of the condensate. 

Comparing the above equations with that of a pure 
harmonic potential, we find the emergence of two new 
terms in Eq.JSJ and ©, i.e., the third and the fourth 
terms in the left-side. Obviously, the third term is di- 
rectly from the distortion of potential. Whereas, the 
fourth one represents the response of coherent wave to 
the distortion of potential manifesting a coupling between 
the motions of center and width. 



III. COLLECTIVE MODES AND FREQUENCY 
SHIFT DUE TO THE ANHARMONIC 
DISTORTION 

When we consider the contribution from the quartic 
distortion, i.e., A ^ in the potential V(x), the nonlinear 
coupling between the oscillations of the center and width 
emerges. The equilibrium points of Eq.© and Eq.© 
correspond to the stable or unstable stationary states of 
the condensate. They satisfy following equations, 



Xo + 2X X 3 o + 3\ Xo w 2 = 0, 



(9) 



(4) 



At a given time t, this function defines a Gaussian dis- 
tribution centered at the position x with width w. The 
other variational parameters t], a and (3 are all real vari- 
ables. Inserting (@J into © we can calculate an grand 
Lagrangian by integrating the Lagrangian density over 
whole coordinate space, L = J_ °° £dx. Then, from the 
Lagrange equations, we obtain the evolution equations 
for all variational parameters. 

The dynamical equations for the center-of-mass and 
width of the condensate is derived as follows 



X + X + 2A X 3 + 3X X w 2 = 0, 



(5) 



(1 + 6\ X 2 o)w + 3Xw 3 = ^ + -^. (10) 

w o 

There is only one stable equilibrium point for A > 0, that 
is, 



Xo = 0, 



wZ wi 



(11) 



(12) 



iu + w + 3Xw + 6X\ iv 



(G) 



where the effective interaction p = —S=, which comes 
from the nonlinear interaction between the particles. 

The other variational parameters can be obtained from 
the center coordinate and the width through the equa- 
tions, 



V^\r,(t)\ 2 w(t) = 1, 



(7) 



For A < there are several equilibrium points, one of 
them is stable and others are unstable. The stable equi- 
librium point also satisfies Eci. (|ll|) and Eci. (|12|) . 

Expanding Eq.(|SJ) and Eq.© around the equilibrium 
points defined by (jll|l and (|12|l and making a routine di- 
agonalizing process, we can obtain following frequencies 
for low-energy excitation modes: 
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which are related to the coupled variation of center and 
width of the condensate 0, B H 0, @- When A = 0, 
u>i corresponds to the dipole oscillation (m = 1) charac- 
terizing the motion of the center-of-mass, and ll>2 is the 
frequency of the variation of the condensate width, it is 
just the low-lying collective mode (m = 0). 

When the potential is not perfectly harmonic, i.e., 
A ^ 0, the contribution from the quartic term will give 
rise to a shift on the frequencies. When A > the fre- 
quency will be blue shifted and when A < the frequency 
will be red shifted. This is true for both single particle 
and BECs. However, it is interesting that for BECs, the 
frequency shift is enhanced dramatically by the atomic 
interaction. The frequencies of two low-energy excita- 
tions as the functions of p for above parameters are plot- 
ted in Fig.l. 
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FIG. 1: Frequencies of two low-energy excitations as functions 
of effective interaction p for different A. 



FIG. 2: The equilibrium width wq as functions of A for dif- 
ferent effective interaction p. 



~" 1.2 




From Eq. l|13|) . we see that the contribution to the fre- 
quency shift comes from the second term in the right- 
hang, i.e., ~ Awq, which is due to the response of co- 
herent wave to the distortion of potential. On the other 
hand, the width wq of the wave function will be broaden 
by the atom interaction. In Fig. 2. we show this ef- 
fete by plotting the dependence of wq on A for different 
interaction parameters. 

From the above discussion, we know, although the an- 
harmonic distortion is very small, the frequency shift 
maybe large due to magnification effects from the atom 
interaction. To demonstrate it, in Fig. 3 we plot the fre- 
quencies of dipole motion of BECs wave-packet for dif- 
ferent atom interactions and anharmonic parameters, ft 
is clearly shown that, the atom interaction will give rise 
to 20% or more shift of frequency when there is only 2% 
of the anharmonic distortion (see Fig. 3, data at p = 20). 
We expect these phenomena can be observed in the fu- 
ture experiment. 



FIG. 3: Frequencies of mass center motions in anharmonic 
trap as functions of A for a classical single particle and BEC 
with different atom interactions p. 



IV. COLLAPSE AND REVIVAL OF THE 
COLLECTIVE EXCITATIONS 

Another promising direction is to investigate the col- 
lapse and revival of the collective excitations [TEl IT^ |. 
which are directly induced by the nonlinear coupling 
effect between two oscillation modes. In the previous 
works, the nonlinear coupling originates from the intrin- 
sic interaction between particles in the system. Here, we 
show a different mechanism for the nonlinear coupling 
that is due to the nontrivial anharmonic corrections to 
the trap. 
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equilibrium point, we have 
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FIG. 4: Variations of the condensate width w for the pa- 
rameter p — 0.4 and the different values of the anharmonic 
parameter A. From up to down A=0, 0.02, 0.04, 0.05, 0.053, 
0.06, 0.08. 



In our case, the qualitative results can be obtained 
with making analysis on Eqs.© and ©. The oscillation 
of the center of the condensate couples with the motion of 
width through the anharmonic parameter A. It is noted 
that the coupled motion can be triggered by putting a 
small shift on the center of the condensate from its equi- 
librium point. The consequent motions of the width of 
the condensate excited by the shift are illustrated in Fig. 
4 for different anharmonicity parameters. In all the above 
calculations, initial shift is set as Ax = 0.1. Fig. 4 clearly 
shows the collapse and revival of oscillation patterns for 
the condensate width with respect to time, induced by 
the anharmonicity. 

Moreover, the change of collapse and revival are not 
monotonic with increasing A and the revival period can 
be effectively controlled by adjusting the anharmonicity 
parameter. In fact, by linearizing Eqs. (JSJ and 10 around 



— QXwqx' 2 - Obviously, the motion of width behaves like 
a periodically driven oscillator. Since the frequency of 
the 'external force' is 2u>i, and the intrinsic frequency of 
the width oscillation is LU2, so the linear combination of 
the two frequencies gives the frequency of the collapse 
and revival, that is, \2uji — Detailed analysis also 
suggest that the frequency of the collapse and revival 
increases monotonically with the nonlinear parameter p, 
but is not monotonic with increasing A. It will vanishes 
for some parameters, e.g. at A ~ 0.053 for p = 0.4. This 
is confirmed by our numerical simulations, as in Fig. 4, 
around A ~ 0.053 for p — 0.4, the period of collapse and 
revival becomes much longer. 

V. CONCLUSION 

We have investigated collective excitations of a Bose- 
Einstein condensate in an anharmonic trap with using 
variational approaches and obtained the analytical ex- 
pressions for the frequencies of the low-energy excita- 
tions. It is shown that the two low-energy excitation 
modes, corresponding to the variations of the center 
and width of the condensate, couple with each other. 
The blue-shift and red-shift on the excitation frequency 
caused by the anharmonic distortion is revealed and 
found to be more dramatic in the case of strong atomic in- 
teraction. Furthermore, the collapse and revival of collec- 
tive excitations in the anharmonic potential is discussed. 
We hope our theoretical results will stimulate the exper- 
iments in the direction. 
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